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The detection of High Energy Cosmic Rays (HECR) with energies around and be- 
yond GZK expected cutoff has introduced the idea of existence of a decaying Ultra 
Heavy Dark Matter (UHDM). If this type of particles make a substantial part of 
the CDM, their decay can have important implications for evolution of the large 
structures and high energy backgrounds. Here we present preliminary results of 
numerical solution of Boltzmann equations in presence of a decaying CDM. The 
evolution of baryons and photons energy are calculated. We find that in a homo- 
geneous universe, UHDM can not have a large contribution to the CDM. However, 
this hypothesis can probably survive if one takes into account the dumpiness of 
the matter in the Universe. 

1 Motivations and Evidences for a Very Heavy Dark Matter 

One of the recently raised and yet unexplained issues in Particle Physics is 
the observation of Ultra High Energy Cosmic Rays (UHECRs)(For a review 
see The source of these particles can not be in a distance larger than 

30 — 50Mpc, otherwise, interaction with CMBR and IR photons considerably 
reduces the energy of the primary particles EH. This phenomena should pro- 
duce a steep cutoff (GZK cutoff) in the spectrum of CRs at energies around 
~ 10 18 eU to l0 2Q eV. Observation shows a local minimum in the spectrum 
around these energies, but unexpectedly it rises again at higher energies. In 
fact the spectrum is compatible with two distinct class of sources, one for en- 
ergies less than ~ 10 18 eU, and the other for higher energies. It is very difficult 
to detect the composition of UHECR primaries. It seems however that the 
data is more compatible with proton or neutrino like primaries Ej. 
No candidate source has been observed in the directpn of UHECRs in a dis- 
tance permitted by GZK cutoff. N. Hayashida et al.h3 reports a small excess 
of ~ 4% for UHECRs with energies ~ 10 18 eU in the direction of the Galactic 
Center and Cygnus region. However, it is shown 113 that there is not any cor- 
relation between observed events and any of individual galaxies in the BOAlpc 
radius for events with energies > 10 18 eU. 

Standard production mechanism of UHECRs is the acceleratation of charged 
particles in the shock waves of AGNs, SNs, or in-falling gas in rich clusters. 
The main condition for these objects to be able to produce the observed spec- 
trum of UHECRs is the existence of relatively strong magnetic field with large 
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coherence length. The secondary acceleration of charged particles by remnant 
of old SNs can accelerate protons up to ~ 10 15 eV and iron to ~ 10 18 eUEZl. 
Models for proton acceleration in AGNs can produce UHECRs up to energies 
~ 10 20 eV c3 with a marginal maximum energy of 4 x 10 21 eUti They also 
predict a change of composition from iron to proton. However, the uncertainty 
on jet density, coherence length of the magnetic field in hot spots and large 
distance to the candidate sources incompatible with energy loss due to inter- 
action with CMB, make these sources marginal candidates. 
Gamma Ray Bursts (GBR) also have been proposed as the candidate source for 
UHECRs E£l. If UHECRs are produced in a relativistic shock wave, one expect 
to observe after-glows of this process in lower energies (like GBRs). Attempts 
to find an excess oi TeV range 7— Rays in the direction of observed UHECR 
events have failed 0. One more argument against an acceleration source for 
UHECRs is that one expects to observe a bunch of UHECRs at least for the 
most nearby sources as it is the case for GRB. No such event has been yet 
observed. 

Another possible class of sources for these events is the decay of a very heavy 
elementary particle with a GUT-scale mass 13. These particles can be either 
long life very heavy particles, or short life particles produced from the decay of 
topological defects. In the latter case, the heavy particles decay in their turn 
to ordinary particles and make the observed UHECRs. There is a number 
of mechanisms which can produce long life ultra heavy particles in the early 
universe. At the end of an inflationary epoch, the parametric resonance of 
inflaton field oscillation .can produce very heavy particles with masses larger 
than inflaton mass itselflia. In super-symmetric models where SUSY breaking 
in the hidden sector is communicated to the observable sector by gauge fields, 
the existence of a global symmetry can lead to stability or long life of lightest 
particle. The lightest messenger gauge field get^lso a large mass and under 
some symmetry conditions can have a long life E3. Vacuum fluctuation of a 
background field coupled to the gravitation also can produce very heavy par- 
ticles with fix ~ 1 for a range of inflaton mass and self couplingp. This makes 
these particles a very good candidates for_Dark Matter (DM). The unitarity 
constraint on the mass of DM particles E3 can be turned over if they have 
never been thermalized. The condition of a long life time for.-a very heavy 
particle can be fulfilled if it has some discrete gauge symmetry EJ. The possi- 
bility that topological defects be the source of UHECRs has been also studied 
extensivelyli3. Some years ago, there was a lot of hopes and attempts to prove 
that topological defects can be the seed of LSS and CMB fluctuations (see H 
for a review). However, the results of comparison between simulations and 
observed power spectrum of LSS and CMB anisotropy practically rules out 
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this possibilityB. 

By contrast, as mentioned above, .some models can produce a CDM totally 
composed of ultra heavy particles Their decay can have important impli- 
cation on the evolution of LSS. Remnants contribute to the total amount of 
baryonic matter in the universe mostly as a non-thermalized hot component. 
Another important implication is the production of large amount of high en- 
ergy backgrounds. This can be used as a test for verification of this hypothesis. 
In this work, we calculate the opacity of the universe for remnants with high- 
est energies and the evolution of their energy spectrum. We compare expected 
flux if CDM is UHDM, with observations. 

2 Model 

We assume that the X particles are neutral bosons with a mass mx = W 24 eV 
and a life time r = 1.2 10 yr. If the decay of these particles looks like the 
decay of Z°, theoretical and experimental arguments show that leptonic and 
hadronic channels have a branching ratio of ~ 1/3 — 2/3. Other channels have 
a negligible contribution. As the dominant channel is the hadronic one and 
the leptonic channel makes very soon a cascade, here we only consider this 
mode of decay. It is very likely that X particles don't decay directly to known 
particles. This leads to softening of the energy spectrum. For simulation of 
hadronization, we consider the hadronization of two gluon jets. Experimental 
observation as well as MLLA (Modified Leading Logarithm Approximation)^, 
LPHD (Local PartonJTadron Duality) ( E3 and references therein) and string 
hadronization model u obtain a softer spectrum with higher multiplicity for 
gluon jets than for quark jets. This is probably more convenient for the decay 
of X particles. 

We use PYTHIA program @ for jet hadronization. This program, like 
many other available ones, can not properly simulate ultra high energy events, 
not only because we don't know the exact physics at l0 16 GeV scale, but also 
because of programming limits. For this reason, we had to extrapolate sim- 
ulation results from lower energies. Fig.0 shows as an example, the photon 
multiplicity for hadronization of a pair of gluon jets. At energies larger than 
l0 12 eV, the multiplicity per log(E) ~ y (rapidity) is roughly constant. The 
same behavior exists for other species. It is the result of KNO scaling predicted 
by MLLA. We use this property to extrapolate the multiplicity spectrum to 
invariant CM mass of 10 24 eV. It should be mentioned that at high energies, 
one expects that KNO scaling be violated. The fact that in Fig.l we don't see 
this effect shows that at high energies, the results of simulation programs can 
be uncertain. We think that various parameterizations of hadronic spectrum 
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Figure 1: Photon multiplicity 

— 1 10 20 30 m hadronization of a pair of 

Log Energy eV gluon jets with CM energies from 

i0 12 eV to 10 24 eV. 

also suffer from the same problem, simply because they have been fitted to low 
energy models or observations. In this respect, simulations are probably more 
close to reality because they take into account more non-perturbative effects. 
In the simulation, we assume that all particles except e , p ± , v, P&7 will decay. 
The mean energy contribution of the stable species is as the following: 13% e^, 
25% p ± , 36% vhv & 25% 7. More than 99% of the total energy belongs to the 
particles with energies higher than 10 20 eU. If baryogenesis happens in GUT 
scale, we expect an excess of baryons over their anti-particles. We assume a 
positive baryonic number and a null leptonic number. This results to an excess 
of e~ and p + over e + and p~ (to keep the universe electrically neutral), and 
v over v. For simplifying our simulation, we assume that all u' s are v e and 
massless. It is interesting to note that z/ s have the largest contribution to the 
total decay energy. 

The decay of a very heavy particle has two implications: Energy dissipation of 
the remnant in one hand heats the ordinary (visible) matter in the Universe 
and on the other hand it contributes itself to the yield of the visible matter. To 
calculate the spectrum of the remnant, one has to solve the complete Einstein- 
Boltzmann equations. Even in the case of linearized equations, the presence of 
interactions makes this system of equations non-linear and coupled. Specially, 
in the case of UHECR energy dissipation, the problem is highly stiff in the sens 
that one has to treat 29 orders of magnitude in energy, from CMB energy to the 
mass of X particles. The numerical solution of these equations is under testo. 
Here we discuss the results of a more simplified model. We consider a homoge- 
neous universe and we neglect the production of high energy particles due to 
interaction of remnants with backgrounds. The reason behind this assumption, 
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Figure 2: Optical depth versus en- 
ergy at different red-shift. Ener- 
gies are particle energy at z = 0. 
Solid z = 0.02, long dashes z = 
0.1, dash 3-dot z = 0.5, dash dot 
2 = 2, dashed z = 10, dot z = 50. 



apart from simplifying the model, is that we are only interested in comparing 
the predicted flux for very energetic particles E > 10 17 eV with observation, as 
we know that lower energies have other sources than presume UHDM and thus, 
they are not directly signals of these particles. We assume that particle energy 
completely dissipates and makes a high energy background. For 7 and neutri- 
nos, we assume that this background has a Gaussian distribution. For massive 
particles i.e e and p ± , we assume that they have a fixed mean temperature. 
In one hand this simplification is motivated by the results of S. Lee iro where 
propagation of UHECRs has been studied for a flat homogeneous universe. He 
considers also the effect of production of high energy particles from interaction 
of remnants with backgrounds and show that there is an increase of photon 
background at energies lower than ~ 10 14 eV. On the other hand, there are 
numerous evidence of existence of a hot IGM, ICM (intra-cluster medium). 




Figure 3: Energy Flux at z = 
with 4 different 7 backgrounds: 
Solid and dotted E bg = 10 10 eV, 
dashed and dash E bg = 10 6 eV, 
and respectively with and without 
radio background. Dots present 
observed flux. Only highest ener- 
gies should be compared with ob- 
servation. 



Backgrounds else than remnants are BBN baryons, primordial neutrinos and 
CMB and star light. We use models developed by Salamon & SteckcrE-3 and 
by Vernet & Valls-GabaudES respectively for IR and visible and for UV back- 
grounds. We assume f2 m = 0.3 and Oa = 0.7. 

As expected, the optical depth depends on the temperature assumed for the 
remnant background. Neutrinos temperature is assumed to be 10 12 eV. The 
reason for this choice is that at lower energies their cross-section is very small. 
Remnant photons have significant role on the opacity of high energy particles. 
We studied following temperatures: 10 4 ey, 10 5 eV, 10 6 eV and 10 10 ey. Each 
of these temperatures has its maximum effect for a range of remnant energies. 
Fig.|| shows as an example, the optical depth for E^ g = 10 6 eV. Reducing 
background energy increases the optical depth because the number of photons 
will increase. 

Fig.|| shows the density of stable particles at z = 0. Irrespective of background 
photon energies, the expected flux of UHECRs at z = from the model is much 
larger than observed flux. This result definitively rule out the existence of large 
contribution of UHDM in the CDM or its life time must be at least 4 order of 
magnitude larger than the age of the Universe (see alsou). However, we believe 



6 











5 



10 
,-10 




10 



-15 



o 8 io 1 9o 1 W% 1 % 2 9o 22 

Log (Energy eV) 



Figure 4: Optical depth of a 
clump with radius = IMpc and 
— = 200 for and 7 solid line, 

dash dot, u dash 3-dot and v 
long dashes. 



that our model is too simplified. The Universe is far from having a uniform 
distribution of matter. Dark Matter as well as baryonic matter is clumpy. The 
opacity for remnant is proportional to the density of the background. Fig]!] 
shows the optical depth of a clump with a mean over- density of ^ = 200 and 
a radius of IMpc. Apart from vs, all other high energy species dissipate their 
energies in the clump. This results to a large reduction of expected flux on the 
Earth and can probably survive the UHDM model. More results about the 
effect of dumpiness on the UHECRs will be reported elsewhere. 
I am grateful of J. Vernet and D. Valls-Gabaud for providing me with the 
results of their UV background model. 
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